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Abstract: The material below the crust of a neutron star is understood to be de- 
scribable in terms of three principal independently moving constituents, identifiable as 
neutrons, protons, and electrons, of which the first two are believed to form mutually 
coupled bosonic condensates. The large scale comportment of such a system will be 
that of a positively charged superconducting superfluid in a negatively charged "nor- 
mal" fluid background. As a contribution to the development of the theory of such a 
system, the present work shows how, subject to neglect of dissipative effects, it is pos- 
sible to set up an elegant category of simplified but fully relativistic three-constituent 
superconducting superfluid models whose purpose is to provide realistic approximations 
for cases in which a strictly conservative treatment is sufficient. A "mesoscopic" model, 
describing the fluid between the vortices, is constructed, as well as a "macroscopic" 
model taking into account the average effect of quantised vortices. 
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1 Introduction 



According to the generally accepted understanding, apart from an outer magnetosphere 
of relatively negligible mass and perhaps also, in the most massive examples, of an inner 
core consisting of matter in some more or less exotic state at more than nuclear density, a 
neutron star can be described in terms of three principle layers: as well as the degenerate 
gas of electrons (and muons) that is present throughout, the outside layer consists of a 
"dry" solid upper crust formed just by a lattice of positively charged ions in a degenerate 
electron gas; the middle layer consists of a "wet" lower crust, whose solid ionic lattice is 
interpenetrated by a neutron superffuid; finally the interior region consists of a neutron 
superfluid in which the ionic lattice has dissolved to form a superconducting proton 
fluid. The purpose of the present article is to develop a simplified model that can 
be used to represent the dynamics of this superconducting superfluid interior, which 
constitutes the greater part of the mass of a typical middlesized neutron star. 

For purposes for which only a very crude description of the bulk motion of the in- 
terior superconducting superfluid region is needed, it may be sufficient to use a model 
involving only two independently moving constituents, of which one represents the 
neutron superfluid while the other represents everything else, since short range elec- 
tromagnetic interactions will tend to ensure that the positively charged protons and 
the negatively charged electrons will tend to move together as an approximately rigid 
background. One of the advantages of such a two constituent description is that it is 
also applicable to the crust layers, so that - in an approximation in which details such 
as the small anisotropies due to the solidity of the crust are neglected - it can be used 
globally to provide a continuously unified description of the star as a whole. A suitable 
two- constituent fluid model for this purpose has recently been constructed |l] in a rel- 
ativistic form that includes allowance for "transfusion", meaning the slow transfer of 
baryonic matter (microscopically attributable to processes such as beta-decay of neu- 
trons to protons) between the neutron superfluid and the other "normal" constituent. 

Although adequate for many purposes, the two-constituent fluid treatment [1]] that 
has just been alluded to is insufficient for the detailed description of phenomena such as 
rotation frequency "glitches" - which are observable via the pulsar phenomena - that 
depend on the small deviations from isotropy that occur not only in the crust (as a 
consequence of its solidity) but also in the interior, due to the effect of vorticity quan- 
tisation and also (at least in the cases that are observable via the pulsar phenomena) 
to the presence of large scale magnetic fields. Using an variational approach based on 
the use of a Kalb-Ramond gauge potential in the manner originally introduced for the 
treatment of ordinary barotropic perfect fluids 0, the present authors have previously 
shown how to construct a kind of model that - at least over short and moderate 
timescales for which "transfusion" processes can be neglected - can provide a satisfac- 
torily relativistic macroscopic description of the effect of the anisotropy due to fibration 
by quantised vortices in what is describable, at a "mesoscopic" (submacroscopic but 
more than microscopic) scale outside the vortices, by a simple superfluid model. How- 
ever although this kind of anisotropic superfluidity model is directly applicable to the 
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terrestrially familiar example of liquid helium, the utility for neutron stars of such a 
single-constituent representation is limited by its neglect of the electromagnetically in- 
teracting constituents, namely the electrons and the protons, which together represent 
only a small fraction of the mass density, but which nevertheless play an essential role 
in phenomena involving magnetic effects. 

It has already been shown that the variational treatment in terms of Kalb-Ramond 
potentials can be straightforwardly generalised from the original single-constituent per- 
fect fluid case p| to the case of perfect two-constituent electrically conducting fluid 
models (including the special kind - with the negatively charged constituent re- 
stricted to be irrotational - that can be used for a "mesoscopic" description of the 
experimentally familiar kind of type II electronic superconductor using a crude but 
fully relativistic representation that neglects the anisotropy due to the solid, not fluid, 
structure of the positively charged ionic background constituent). Following this exam- 
ple, the purpose of the present work is to show how - still neglecting the "transfusion" 
processes that are likely to be important only in the very long run - our previous macro- 
scopic model||, representing the effect of vortex fibration in a simple superffuid, can 
be extended to a three-constituent model providing a description that is applicable to a 
superconducting superffuid in which the other two independently moving constituents 
are the superfluid protons and the degenerate electron gas. 

A considerable amount of work on three-constituent models for this purpose has 
already been carried out by several authors^, ^|, 0, ||, but until now only within the 
framework of a "non-relativistic" Newtonian description (as a generalisation of the kind 
of two-constituent model used for treatment of an ordinary "type II" metallic supercon- 
ductors of the sort familiar in a laboratory context, for which the two constituents are 
superconducting electrons and a solid ionic background lattice). Some of this previous 
work includes allowance for details such as dissipative effects e.g. to various kinds of 
resistivity, that in the present description will for the sake of clarity will be provisionally 
left aside for treatment in a subsequent treatment. The omission of such dissipative 
effects is however by no means the only reason why the conservative treatment given 
here turns out to be technically simpler and more elegant than its "non-relativistic" 
predecessors: the other reason is that the Galilei invariance group applicable to the 
Newtonian description is more complicated than (since it is a degenerate limit of) the 
(in the technical sense "semi-simple") Lorentz invariance group that is locally applica- 
ble to the fully (general as well as special) relativistic treatment provided here. When 
the electromagnetic coupling is taken into account, the relativistic description is even 
more satisfying in that the material part and the electromagnetic part are treated on 
the same footing instead of the awkward mixture of Galilean and Lorentzian invariances 
in the Newtonian context. 

It is to be remarked that the improvement not just in physical accuracy but also 
in mathematical elegance that is obtainable by going from a Newtonian to a relativis- 
tic treatment has already been demonstrated in the case of Landau's original two- 
constituent superfluid model (in which the second constituent just represents the en- 
tropy current), whose relativistic generalisation has provided insights that have lead 
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to a more efficient "canonical" description [Kj of its original Newtonian version. This 
suggests, as a challenge for future work, that the traditional presentation of the more 
elaborate three-constituent fluid models developed in previous work|| [7|, §], should also 
be susceptible to similar improvements in elegance by an analogous "canonical" refor- 
mulation along lines that would presumably be obtainable by taking an appropriate 
Newtonian limit of the even more elegant relativistic kind of model developed here. 

Of greater importance than the analysis of the Newtonian limit, as a challenge for 
future work, the kind of model developed here will need to be completed by the provision 
(on the basis of appropriate microphysical analysis) of explicit forms and values for the 
functions and parameters involved , while it will also be necessary to develop the more 
extended treatment needed to include suitable allowance for dissipative effects (not to 
mention refinements such as allowance for the spinor as opposed to scalar nature of the 
Cooper pairs involved in the microscopic description |TT| of both the neutron superfluid 
and the superconducting proton current). 

Another challenge is the development of a corresponding relativistic description for 
the solid crust layers. An appropriate relativistic formalism for describing the purely 
solid "dry" upper crust layer of a Newton star was developed quite a long time agofT2|, 
[13], |IH , though its technical complexity has delayed its effective application to specific 
astrophysical problems [|15]. What has yet to be constructed is the even more elaborate 
generalisation required for treating the "wet" lower layer of the solid crust: machinery 
capable of providing a treatment on a "mesoscopic" scale of the interaction between the 
permeating neutron superfluid and the solid background lattice has already been made 
available [0] , but it still remains to develop the corresponding macroscopic analysis 
required for the relativistic description of the large scale effects of phenomena such as 
vortex pinning ]T7| in such a "wet" solid. 

As declared above, the ultimate task of the present work is to provide a kind of 
model that can account for the macroscopic anisotropies attributable to quantised vor- 
tex lines associated with both rotation and magnetic fields, in the deep superconducting 
superfluid interior of a neutron star. However, as a preliminary step, we must first de- 
velop a suitably relativistic "mesoscopic" (submacroscopic but not microscopic) scale 
treatment applicable to the superconducting superfluid in between the relevant string- 
like vortex defects. The most elegant versions of the models presented below are those 
that apply to the zero temperature limit, which should be a good enough approxima- 
tion for most relevant astrophysical applications in neutron stars. We shall however 
include a rudimentary allowance for the effect of a finite temperature associated with 
a non- vanishing conserved entropy distribution comoving with the degenerate electrons 
in the negatively charged "normal" constituent. 



2 Three constituent perfectly conducting fluids. 

Following what has long been established as a routine procedure |T8f - (the natural rel- 
ativistic analogue of the phenomenological approach developed in a Newtonian context 
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19| as a generalisation of Landau's original two-constituent model) - the task of the 
present section is to set up the category of simple three-constituent perfectly conducting 
fluid models that is needed as a preliminary basis for the more specific developments 
that follow. As the specialisation to be treated in the next section, this category in- 
cludes the particular restricted case that is appropriate for the treatment of the kind 
of superconducting superfluid that is relevant in neutron stars on a "mesoscopic" scale, 
meaning a scale large compared with that of the underlying microscopic particle de- 
scription, but small compared with the macroscopic scale of separation between the 
vortex defects on which the superfluid comportment is locally violated. On the other 
hand the category set up in the present section is itself a just specially simple limit 
within the more general category to be developed in the later sections for the purpose 
of treating the superconducting superfluid on a "macroscopic" scale, meaning a scale 
that is large compared with the separation between vortices. 

The three independent constituents under consideration are identifiable as the su- 
perfluid neutrons - which make up most of the mass density in the relevant neutron 
star layers - with number current four- vector n p , the superconducting protons, - which 
make up a small but significant part of the mass density - with number current four- 
vector n p , and finally the degenerate non-superconducting background of electrons with 
number current vector n p . The latter make up a negligibly small fraction of the mass 
density, but they nevertheless have a crucially important role in so far as electomag- 
netic effects are concerned, since in terms of the electron charge coupling constant e the 
corresponding total electric current vector four vector will be given by 

f = e(n p - n p ) . (1) 

As well as the three principal constituents that have just been listed, our discussion 
in this section will also include allowance for a fourth constituent, namely the entropy. 
However this will not be treated as fully independent, since (as will be a realistic 
approximation for most of the relevant astrophysical applications) its current four- 
vector s p will be postulated to be constrained so as to be aligned with that of the 
electrons. This means that, apart from the "super" constituents with currents n p and 
n p , there will be just a single independently moving "normal" constituent with current 
in the direction specified by a unit vector u p , as characterised by the normalisation 
condition 

u p u p = -1 , (2) 

(on the usual understanding that the units are such that the speed of light is unity) 
with 

< = n e u p , (3) 

s p = su p , (4) 

where the scalars s and n c are respectively the entropy density and the electron number 
density in the "normal" rest frame specified by u p . 

Since our present treatment will be restricted to the conservative limit in which 
dissipative effects are neglected, the analysis will be assumed to be expressible in terms 
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of a variational principle based on a Lagrangian density that, in the absence of elec- 
tromagnetic effects, can be presumed to be given by a function Am say depending just 
on the six independent scalars obtainable by mutual contractions of the three inde- 
pendent current vectors n£, ng, n£ and on the entropy density s. The independent 
variations of these quantities determine the effective momentum covectors //^ and piP 
respectively associated with the superfluid neutrons and the superconducting protons, 
together with the "normal" momentum covector fi e p associated with the electrons, and 
the corresponding temperature B, via a prescription of the form 

<9A 

<5A M = /i n p 8n£ + fi p p 5nP + // p <5< - 6 5s + — ^ 5 9pa , (5) 



in which, as an elementary Noether type identity, we shall automatically have 

\(jf v n> + n* u n> + pL%n!)!r- (6) 



<9A M = dA M 
dg P a dg^p 



It is useful to allow here for the possibility of varying the background spacetime metric 
g pa , not only for the purpose of dealing with cases in which one may be concerned 
with General Relativistic gravitational coupling, but even for dealing with cases in 
which one is concerned only with a flat Minkowski background, since, as will be made 
explicit below, the effect of virtual variations with respect to the relevant curved or flat 
background metric can be used for evaluating the relevant "geometric" stress energy 
momentum density tensor T pa . 

The preceding formulae (|5|) and (||) can conveniently be combined in a single ex- 
pression of the more concise form 

5 A M = [i% <K -988+ \n\ < g va 8g pa , (7) 

using the summation convention for "chemical" indices represented by capital Latin 
letters running over the three relevant values, namely x=n,p,e. Using this convention, 
the equation (P for the electric current density can be rewritten in the concise form 

f = e x n£ . ( 8 ) 

where the charges per neutron, proton, and electron are given respectively by e 11 = 0, 
e p = e, and e e = — e. 

The standard minimal prescription for inclusion of electromagnetic interactions is 
to use a combined Lagrangian scalar density in which the "matter" contribution Am is 
augmented by an electromagnetic field contribution A F and a gauge dependent coupling 
term of the usual form to give a total Lagrangian scalar £ expressible as 

C = A + j"A p , (9) 

where the gauge independent part has the form 

A = A M + A F , (10) 
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and where A p is the electromagnetic gauge form, while the gauge independent contri- 
bution Ap is given in terms of the corresponding electromagnetic field tensor, 



F pa = 2V [p A a] , (11) 

(using square brackets to indicate index antisymmetrisation) by the standard Maxwellian 
formula 

A F = —F pa F ap . (12) 
16tt p y 1 

The variation of the total will then be given by 
5C = 7T- 5n p -B5s + f 5A p + i-F CT " 5F P<7 + l - (y* < <T + ^ P F V °) 5g pa , (13) 

where the gauge dependent total momentum covectors are given by 

Tr* p = fjX p + exA p . (14) 

In order to characterise a variation principle that can be used to specify the equations 
of motion of a multiconstituent medium, the electromagnetic potential covector A p is 
of course allowed to vary freely, but the variations of the independent currents n p 
cannot be allowed to be arbitrary, since this would evidently lead to an overdetermined 
system in which the momenta would simply have to vanish. The standard procedure for 
restricting the allowed current variations so as to obtain a system with the appropriate 
number of degrees of freedom for an ordinary fluid mixture is to require that the the 
current variations be determined by free displacements of the corresponding worldlines 
- in the manner originally introduced for the particular case of a relativistic perfect 



fluid by Taub pOfl . The ensuing variational equations then automatically ensure the 
conservation under transport by each current flow n p of the corresponding generalised 
vorticity tensor as defined by 



^ pa 



2V [p 7r\] = 2V [p /i x ff] + e*F pa . (15) 



Such an exterior derivative will of course automatically be "closed", i.e. its own exterior 
derivative will vanish: 

VpV]-0. (16) 

Within such a perfect multiconstituent fluid category, the special properties of su- 
perfluidity and superconductivity are characterised by the constraint that the corre- 
sponding conserved vorticity should vanish initially, and hence, since it is conserved, 
throughout the ensuing evolution of the system. Such constraints - which in our present 
application concern the neutrons and the protons but not the electrons - are inter- 
pretable as integrability conditions for corresponding potential scalars. For the neutrons 
we shall have 



w» pa = =► vr n p = ^ p = 5 V p ^ , (17) 
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while for the protons we shall have 

w% = =► ti* = ft + eA p = | V p V p , (18) 

where the locally defined potentials </? n and f p are interpretable as phase angles (with 
period 2n) associated with underlying boson condensates, and the factors 2 in the 
denominators are included to allow for the fact that the relevant bosons are presumed 
to consist not of single protons or neutrons but of Cooper type pairs. 

The way the Taub type "connective" variational procedure works out is as follows. 
For any current vector n p the corresponding variation rule is interpretable as meaning 
that the ensuing variation of its dual three-form, as given - in terms of the antisymmetric 
measure tensor e^ upa associated with the 4-dimensional spacetime metric g pv - by 

should be specified by its Lie derivative with respect to the relevant worldline displace- 
ment vector field, £ p say, i.e. 

<*Vp = Vp + 3K|>«,V p] £ ff , (20) 

where V M denotes the covariant differentiation operator specified by the relevant (flat 
or curved) spacetime background metric g pv . 

When applied to the particular case of the "normal" (i.e. non-superfluid) electron 
current n e p , for a corresponding displacement vector field £ e p , the foregoing prescription 
can be seen to give a variation of the form 

K = cv. < - <v CT a + <v CT c - i < ST 6g,u, , (21) 

from which it follows that the variation of the scalar electron number density introduced 
in (0) will be given by 

5n e = CV ff n e + n e u p u a V a g + n e V a C - \n e {g^ + u»u v ) 5 9flv , (22) 

Since we are postulating that the electron current is comoving with the entropy cur- 
rent as part of a combined "normal" constituent, for which the shared unit 4- vector 
appearing in (|J) and @ will thus be subject to the variation law 

5u p = C CT V CT u" - u°V a £«f - u p u v u°V a C + iu'W Sg^ , (23) 

it follows that the variation of the entropy current will be governed by the same dis- 
placement vector field £ e p as that of the electron current, and hence the variation of the 
entropy density introduced in (|j) will be governed by a rule of the same form as (P%fl, 
namely 

5s = £TV CT s + su p u°V„ £f + sV CT C - 5*GT + "V) 5g pv . (24) 



S 



When metric variations are involved, the application of the variation principle re- 
quires allowance not just for the variation of the Lagrangian scalar C, but also for the 
variation of the relevant spacetime measure, which is proportional to the square root of 
the modulus of the metric determinant \g\. This means that the complete variational 
integrand that one needs to evaluate is the "diamond" variation as defined by 

OC = \\g\\- 1 ' 2 5(\\g\\ 1/2 £) =5C + \Cg» v bg pv . (25) 



By taking the variations of the neutron and proton currents to be given in the same 
way as that of the electron current by corresponding displacement vector fields £ p and 
£jf , the required "diamond" variation of the Lagrangian scalar density given by @ can 
be obtained in the standard form 

0£ = f x P £ X P + (f - ^ V ff F pa ) 8A P + 5g pv + V a K\ (26) 

in which the coefficients /* of the independent displacement displacement vector fields 
are respectively interpretable as the force densities acting on the corresponding 
constituents, while the (symmetric) tensor T pv evidently represents the ordinary ( "geo- 
metric" ) stress momentum energy density. It can be seen that the force densities acting 
on the superfluid neutrons and the superconducting protons will be given respectively 
by 

/; = «p + ^V,<, (27) 

and 

P p = n;wl p + 7r p p V ff n; . (28) 

where w n a and w p ap are the corresponding vorticity two- forms as defined by fllS]). The 
force density acting on the "normal" constituent - consisting of the electrons and the 
comoving entropy current - can be seen to be given by an expression of the not quite 
so simple form 

f e p = 2 U CT V [CT II P] + n p V CT u° + % n c - 0V P s , (29) 

in which the one-form Ii p is the effective momentum density of the normal constituent 
(electrons and entropy) treated as a whole, as given by 

n p = n e vr e p + sQu p . (30) 

The stress momentum energy density tensor can be seen to have the form 

T pa = T M pa + T/ a , (31) 

where the part derived just from the material Lagrangian density contribution Am is 
given by 

?m a = <rt + s®u p u a + ^ M g p a , *m=Am- n>* + sQ , (32) 
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while the electromagnetic contribution has the usual Maxwellian form 

T/„ = ^ (F^F va - \f^F^) . (33) 

Although it is irrelevant as far as the variation principle is concerned, it is to be noted 
for the sake of completeness that the remainder current in the final divergence term of 
(26) will be given by 

11° = 27r^iX ] + 2Qu u ^s^ + ^F™ 5A U . (34) 



It is evident from fl26|) that the requirement of invariance of the action with respect 
to free variations of the gauge form A p leads to electromagnetic source equations of the 
usual Maxwellian form 

V CT F pa = Anj p . (35) 

Subject to this condition, it follows as a Noether type identity (obtainable by identifying 
all the displacement vectors with a single vector field with respect to which the system 
is displaced as a whole ITB, UM) that, independently of the other field equations one 



must have a force balance relation of the form 

n + n + n = v, v p , m 

in which the right hand side represents the total force density acting on the medium, 
which will of course vanish for the kind of isolated strictly conservative model with 
which we are concerned here. More particularly, it is evident from (p6|) that the con- 
tent of the Taub type "connective" variation principle expressed by the condition that 
independent adjustments of the displacement fields should have no effect is equivalent 
to the dynamical requirement that each of the separate force density contributions on 
the left of fl36|) should vanish separately: in short the equations of motion of the system 
will be expressible simply by the condition 

n = (37) 

for each of the three relevant chemical index values x=n,p,e. 

Assuming that (unlike what would be required for a model of "transfusive" typejlj) 
there is no algebraical constraint on the momenta, it can be seen (by contracting the 



relevant force formulae (27) and (28) with the corresponding current vectors) that for 
the cases of the neutrons and the protons the conditions ( |3"7D include the implication 
that their currents are separately conserved, 

V p < = 0, V p n£ = 0. (38) 

The remaining contents of the neutron and proton equations of motion will simply 
reduce to the respective forms 

Kwl p = , n>P p = , (39) 
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where the vorticity two forms are as defined by (|l"5|). It is well known that the form 
of such equations ensures that, in view of the closure property flTfi|), the neutron and 
proton vorticity forms will be conserved in the sense of being dragged along by the 
corresponding flows [M], which means that they are consistent with the special require- 



ments for superfluidity and superconductivity respectively, which are that they should 
simply vanish, in accordance with ( |17D and ([18]). 

For the non-superfluid "normal" constituent, the implications of the zero force con- 
dition (|37|) are not so simple. In this case, it can be seen from ( f29"|) that contraction 
with the relevant unit flow tangent vector will just give 

tAr c p V CT < = 9V^ ff . (40) 

To obtain separate current conservation laws that are required, it should be understood 
that the electron current is not specified freely but is given by a procedure of the kind 
appropriate for non-barotropic perfect fluid theory f22fl via a projection (of the kind that 
is standard as the basis of relativistic elasticity theory |12|, p"3[| ) onto a three dimensional 
base manifold endowed with its own measure. This means that in terms of the space- 
time pull-backs of suitably chosen base coordinates X A , A=l, 2, 3, it will be given, using 
the abbreviation X A = V p X A , by an expression of the form n e M = e^ upa X 1 V X 2 p X A a . 
The variation law fl2"T| ) is equivalent to what is obtained by allowing free variations 
of the X A considered as a set of three independent scalar fields on spacetime: such 
variations are equivalent to the specification of a corresponding displacement vector 
field £ e p for which one will have 5X A = £ C P V P X P . Such a construction automatically 
ensures that electron current will be identically conserved, and the dynamical equation 
(p0|) then entails a similar condition for the entropy, i.e. one ends up with the separate 
conservation laws 

V CT < = 0, V ff s" = 0. (41) 

Under these conditions, the remaining contents of the equations of motion of the "nor- 
mal" constituent, as expressed by the vanishing according to ([37]) of the force density 
f e p given by fl2"§|), will take the form 

« p + s CT V ff (Qu p ) + sV p 6 = , (42) 

which - unlike the corresponding dynamical equations (^) for the neutrons and the 
protons - does not have the form of a generalised vorticity conservation law, except in 
the zero temperature limit, G = 0, for which only the first term remains. 



3 Mesoscopic superconducting superfluid model. 

Within the category of perfectly conducting fluid motions, as characterised by the 
equations of motion (0), (|39"D, ([41]), (|42|), the "mesoscopic" superconducting superfluid 
case relevant on an intermediate scale (in between the vortices) will be characterised 
by the constraints ( |TTD and (|i~8"|), whose compatibility with the dynamical equations is 
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ensured by the particular form of flog) . Instead of obtaining them as a particular case 
within the general perfectly conducting fluid category, we can proceed directly to the 
specialised superconducting superfluid equations of motion, as given just by fl3~8|) , (f4TD, 
(^) together - in accordance with ([H]) and (|T8p - with 

2n n p = hV p ^, 2vrP = W p ^ p , (43) 

by using a variational procedure of a modified kind. Since this new kind of procedure will 
treat the "super" currents on a rather different footing from the others, its presentation 
will be facilitated by the introduction of a correspondingly restricted chemical index, 
for which we shall use the capital Greek letters T, $ which differ from the unrestricted 
chemical indices X, Y in having a range that is restricted to exclude the "normal" 
current values. For the three-constituent category under consideration here, the only 
"normal" value to be excluded is X=e. It follows that, for this category, the Greek 
letters will be interpretable as isotopic indices, taking the pair of values T=n for the 
neutron case and T=p for the proton case. With this convention, the two separate 
equations ([43]) can be combined into the more concise form 

27r- = W p ^ r , (44) 

in which it is to be recalled that the factor 2 is present in order that the left hand side 
should represent the effective momentum covector of the relevant bosonic unit, which 
will consist not of a single baryon but of a Cooper type pair. 

According to the rules of the modified procedure, instead of being restricted to satisfy 
displacement variation rules of the kind fl21~|) satisfied by the "normal" constituent, it is 
permissible for the isotopic doublet formed by the "super" currents and to vary 
freely apart from the restriction that the corresponding conservation laws (|3^) should 
be satisfied. The enforcement of these conservation laws, which can be combined in the 
single expression 

V ff n; = 0, (45) 

is to be obtained in the usual way by the introduction of corresponding Lagrange 
multipliers. Provided they are introduced with the appropriate normalisation, these 
multipliers turn out to be identifiable with the potentials </? n and tp p introduced in (|T] 
and (HH). This is done my making the replacement 



C » A , (46) 
where the modified Lagrangian scalar is given by the expression 

C 1 = C + \(p r V„n°, (47) 

in which summation over the isotopic index T is of course to be understood. The 



correspondingly modified replacement of the "diamond" variation formula ( ^6|) is 
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OA = (tt* - |V P ^)5< + / e p + (/ - ^ V CT F") M p 

+ |(V ff <)<V r + iT^ + V CT 7^ , (48) 

in which the force density f* acting on the normal constituent is given by exactly the 
same formula (p9|) as before. Although it is irrelevant for the purposes of application 
of the action principle, it may be noted that instead of fl3~4]), the remainder current in 
the final divergence term of (|48|) will be given by the rather different formula 

Tl{ = \ip r {8nZ + l -n;cr5g^) + 2U u ^ a u^ + ^F™ 8A U . (49) 

More to the point for practical purposes, the stress energy momentum density tensor 
T^ v in (f4"8| ) will be given by the same formulae (|3T|), fl32|), ( p3|) as before, apart from 
the fact that in the new version the relevant generalised pressure function will be given 
by 

*m = A m - ~<V ff <p n + n°(eA a - |V CT </) - </4 + sQ , (50) 

which can be seen to be equivalent to the original version, as given by (|32|), when 
the superconducting superfluid field equations (|43|) are satisfied, as is necessary in this 
modified formulation for invariance of the action with respect to arbitrary infinitesimal 
variations of n£ and n?. 

The essential difference between a "normal" and a "super" current is not just that 
the latter is irrotational so that its momentum covector is locally proportional to the 
gradient of a scalar field, but more particularly that this scalar field will be globally 
definable as a periodic phase angle, so that the corresponding circulation integral will be 
quantised. Specifically, according to ([0]), the relevant neutron and proton circulation 
integrals for a given circuit will be given in the present case by 

j^dx" = k?iv y , (51) 

for a pair of integral proton and neutron phase winding numbers respectively given by 
v n and v p . These numbers will remain constant as the circuit is continuously displaced, 
except when it crosses what, on a mesoscopic scale will be describable as a string like 
singularity, but what on a microscopic scale will be a vortex type defect whose core 
consists of a region in which the superfluidity property must break down. It is expected 
that, in most of the relevant neutron star layers, conditions will be such as to ensure 
that such vortex defects will generally be of elementary type, meaning that they will 
either be simple neutron flow vortices, as characterised, for an appropriate choice of 
orientation, by winding numbers v n = 1, v p = 0, or simple proton vortices, as similarly 
characterised by v n = 0, v v = 1. In order to obtain a description that is applicable on 
a macroscopic scale, it will evidently be necessary to average over a large number of 
such (elementary, or in more exotic circumstances, higher order) vortices. Following the 
example of our treatment of the single constituent case [[J, the purpose of the next two 
sections is to develop and show how to apply the kind of formalism needed to achieve 
this. 
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4 Macroscopic allowance for vorticity fibration 



The standard (Taub type) flow line variational procedure illustrated by the category of 
models set up in Section |2| has the advantage of being extremely versatile: for example 
it can be generalised for application to simple elastic solids 1131) as well as to conducting 



solids | IE||, and it has been used most recently for setting up the "transfusive" two- 
constituent model |] that has been designed to provide a crude but globally applicable 
description of a neutron star as a whole. However for the purpose of dealing with the 
macroscopic effects of vorticity quantisation, we have found that it is more convenient 
to start off on the basis of a very different kind of variational procedure. 

The standard procedure, as used in Section |2|, for the treatment of any current, 
n cr say, is based on a projection of the four dimensional space-time background onto a 
three-dimensional material base space, whose local coordinates can then be pulled to 
specify corresponding comoving space coordinates on the space-time background: the 
dual three-form, *n upa , is obtained as the pullback of a prescribed volume measure on 
the base space. In the alternative variation procedure, which has been shown to lead 
ultimately to the same equations of motion, both for a simple perfect fluid [0] and for 
a perfectly conducting two-constituent model 0], an analogous projection construction 
is used to specify not the current three-form *n upa but the corresponding vorticity two 
form w pa , which is to be obtained as the pullback of a prescribed area measure on a 
two-dimensional base space. In terms of suitably chosen local coordinates x 1 ; X 2 on 
the vorticity base manifold, the corresponding pair of scalar fields induced on the four 
dimensional spacetime background will specify the vorticity according to the formula 
w pa = 2x \ p X ! a ]- Such a prescription automatically ensures that the vorticity two-form 
will satisfy both the algebraic degeneracy condition 

W[^w pa] = , (52) 

and the closure condition 

V[„ w pa] = . (53) 

For the purposes of the variational principle it is to be postulated that the projection 
of the four dimensional space-time background onto the two dimensional vorticity base 
space should be freely variable, which is equivalent to the postulate that x 1 an d X 2 
should be considered as freely variable local fields on spacetime. As for the three 
dimensional projection used for the direct specification of a current, so also for the two 
dimensional projection used for the specification of a vorticity form, it is convenient 
to represent an infinitesimal variation via an infinitesimal displacement vector field, £ p 
say, in terms of which the corresponding base coordinate variations will be given by 
fix 1 — i p X\ P i fix 2 — i p X 2 P - As in the case of currents, the use of such displacement 
vector fields to specify the corresponding "convective" variations enables us to employ 
a formalism that is manifestly invariant with respect to changes of the base coordinates, 
whose explicit use can thereby be avoided in the work that follows. It can easily be seen 
that a base coordinate independent expression for the infinitesimal variation induced 
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by the displacement field £ p will be given by the formula 

5w^ = -2V { ^w u]p e) ■ (54) 

This is evidently interpretable, taking account of the closure property fl52|), as the 
condition that the infinitesimal variation of w pv should be given simply by its Lie 
derivative with respect to the 

The foregoing procedure makes it possible to treat vorticity as an independent field 
in its own right. However in such an approach the essential identity of the vorticity 
two form as the exterior derivative of a corresponding momentum one-form must be 
ultimately be established as a dynamical equation, instead of being imposed in ad- 
vance as a defining property in the manner used in the preceding sections. In order to 
obtain the required exterior differential relation between momentum and vorticity in 
a variational treatment, it has been found convenient 0, |], [5] to replace the standard 
three-dimensional base projection procedure used (as described above) for the specifi- 
cation of the relevant current, n p say, by an alternative Kalb-Ramond type procedure 
whereby the current is specified in terms of a freely variable gauge two-form, i.e. an an- 
tisymmetric tensor field B pv = —B Ufl in terms of which the relevant current three-form 
is given by 

*?Vp = 3V[ M B vp ] . (55) 

By the Poincare lemma and its inverse, the existence of such a gauge potential is a 
sufficient and locally necessary condition for the closure property 

V] = > (56) 
which is equivalent to the conservation law 

V p n" = 0. (57) 

The equivalent statement in dual terminology is that this conservation law fl5"T|) is the 
necessary and sufficient condition for the local existence of a bivector field b pu = —b UfJi 
such that 

n" = V CT b pcJ , (58) 

where b pu is the dual of a corresponding Kalb-Ramond form in terms of which it will 
be given by 

b^ = \e^B pa . (59) 

It is to be noted that a Kalb-Ramond representation of the form (|55|), or equivalently 
(|58|), can not be used in cases such that of the globally applicable "transfusive" model 
in which the relevant currents are not separately conserved. It is also to be remarked 
that (assuming it is admissible) such a representation will evidently not be unique: it 
is subject to gauge transformations specifiable in terms of an arbitrary one-form a v by 

>-> + 2V [fl a v] . (60) 
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In the particular case of the three constituent models with which the present discus- 
sion is concerned, the considerations presented above lead us to introduce an isotopic 
doublet of free gauge bivector fields b? a and h^ u in terms of which the neutron and 
proton currents will be given by the combined expression 

< = V CT £>f, (61) 

where, as in the previous section, it is to be understood that the isotopic index runs 
over the values T=n and T = p. However since the problem of vorticity quantisation 
does not arise for the normal constituent we shall continue to treat it as before, meaning 
that n e p will be supposed to be defined, via a projection, as the dual of the pullback of 
the measure on a three dimensional base space. Treating the gauge fields and b? a 
as freely variable, but treating the associated vorticities w n and w p pa , via independent 
projections, as pullbacks of the measure on a two dimensional base space in the manner 
described above, in a new Lagrangian obtained from (§) by a replacement 

£ ^Cjl, (62) 

where (summing as usual over the isotopic index T) 

C I = C + fe'w*, , (63) 

it can be verified that one obtains a variational formulation leading back to field equa- 
tions precisely equivalent to those obtained by the original procedure in Section |2|. 

The advantage of the new procedure is that, following the example of our single- 
constituent prototype 0, it can easily be generalised so as to provide a macroscopic 
description of the effects (whose analysis in a Newtonian framework was pioneered 
by Bekarevich and Khalatnikov[|2"3"|, |23j) to be expected as a consequence of vortex 
quantisation, by replacing the original formula (pTOj ) for the gauge independent action 
contribution A in the new Lagrangian 

C^A + fAp+'-b^w^, (64) 

by a more general formula according to which A is given as an algebraic function not 
just of the (gauge independent) fields ng, n p , ng, s, F pu , and of course the metric g^ u , 
as before, but also on the (similarly gauge independent) doublet of vorticity two-forms, 
w n prT and w p pcr . This means that its most general infinitesimal variation will be given by 
an expression of the form 

1 A 

5A = f p 5n p + Sw^ + ^ c p Sn p - Q5s + —H ap 5F pa + — 5g pa , (65) 

where, as in (|5|), the coefficients of the metric variations are not independent of the 
others but by the relevant generalisation of the Noether identity (^) must satisfy 



+ \\ v T p wl° + —H vp F v ° . (66) 
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The new doublet of bivectorial coefficients \° p = — X pa in this expansion characterises 
the macroscopic anisotropy arising respectively from the concentration of energy and 
tension in mesoscopic vortices of the neutron and proton superffuids as a consequence of 
their vorticity quantisation conditions, in the manner discussed in our previous work|3| 
on the single constituent model. These new four dimensional bivectorial coefficients 
replace the three dimensional (space) vectorial coefficients introduced for a similar pur- 
pose in a more restricted Newtonian framework by Bekarevich and Khalatnikov||23|| . 
As before, the covectorial coefficients yP p) and \f will be interpretable respectively 
as the mean momentum covectors per particle of the neutrons, the protons, and the 
electrons, and will be interpretable as the temperature of the "normal constituent" . 
Finally the bivectorial coefficient 7i pu = —7i ap will be interpretable as an electromag- 



netic displacement tensor, in terms of which the total electromagnetic field tensor (|TT 
will be given by an expression of the form 

in which M. p(J is what can be interpreted as the magnetic polarisation tensor. In the 
application that we are considering, the - typically dominant - polarisation contribution 
4nAi pa is to be thought of as representing the part of the magnetic field confined 
in the vortices, while the - typically much smaller - remainder T-i pa represents the 
average contribution from the external field in between the vortices, which can be 
expected to vanish by the "Meissner effect" in strictly static configurations, but which 
can be expected to acquire a non zero value due to the "London effect" in rotating 
configurations. 

The analogue of fl26|), as the macroscopic replacement of ((48|), for an infinitesimal 
variation according to the foregoing rules - with the neutron and proton vorticity con- 
gruences convected by displacement vector fields £ p and £A and with the flow lines of 
the "normal" constituent convected by a displacement vector field £ e p as before - will be 
a "diamond" variation (i.e. one that includes allowance for the effect on the background 
spacetime measure) expressible in the standard form 

OCj = (V, K - \wT pa )6b p ° + /* + (f - V ff W°) 5A p + bg pv + V CT K{ , (68) 

in which the force densities /* appearing as coefficients of the three independent dis- 
placement displacement vector fields £ p will be given for the isotopic doublet values 
X=n and X=p by 

r P = « + v, \rx P , r P = k + ^ K v x P . ( 69 ) 

while the third force density f e acting on the "normal" constituent will still be given 
by the same formula Q2§| ) as before. 

Although it is of no relevance for the application of the variation principle, it can 
be noted for the record that the current appearing in the final divergence term of QB"8"D 
will be given by 

K/ = Sb? - (6f + A£>*,g + -ttXV" 8g, v + 2U V ^V ] + 5A V . (70) 
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An entity of much greater practical interest is the corresponding stress momentum 
energy density tensor, which can be seen to be given by 

TJ> = n>* + fl e«X + K"<« + K"<° + + *9 P „ , (71) 

where the new generalised pressure function is given by 

viz = A - n>* + sQ + (V, 7T- - | w y (72) 

The last term in (fT2j) will evidently drop out when we impose the condition of invariance 
with respect to infinitesimal variations of the bivectorial gauge potentials b pcr and b pa 
is imposed, a requirement which can be seen from (|68|) to give field equations of the 
form 

= 2V [p 7r n a] = 2V [P /A] , = 2V \p^\ = 2 Mp/^] + eFpa • (73) 

These expressions are evidently equivalent to what in the previous formulation were 
merely the definitions of the vorticities as given by the general specification (fL5f). 

The remaining field equations obtained from fl6"8|) will consist of an electromagnetic 
source equation having the form 

V a H p<J = 4nj p , (74) 

together with the condition that the force density coefficients should all vanish. For the 
normal constituent, the corresponding condition 

f e p = (75) 
will take the same form (^) as before, but it can be seen from (|69| ) that the other two 



force balance conditions, 

r p = , (76) 
will involve new terms. In particular, the neutron force balance condition can be seen 



from (|69f) to be expressible as 

2<V [CT ^ ]+ < p V,Ar = 0, (77) 

in which the first term is interpretable as the negative of the Joukowski force density 
due to the "Magnus effect" acting on the neutron vortices, while the last term (which 
was absent in the mesoscopic description) represents the extra force density on the fluid 
due to the effect of the tension of the vortices. The proton force balance equation can 
be similarly be seen to be expressible in the form 

In; V [CT fj? p] + en;F ap + < p V, A™ = , (78) 

in which the first term is the negative of the Joukowski force density that acts on the 
proton vortices, the middle term is the Lorentz force density representing the effect of 
the magnetic field on the passing protons, and the final term again represents the force 
density due to the effect of the tension of the vortices. 
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5 Entrainment in partially separated models. 



In order to be utilisable for particular applications, the characterisation of the meso- 
scopic and macroscopic models described in the preceding sections needs to be com- 
pleted by the designation of an appropriate "equation of state" giving the specific form 
of the mesoscopic Am and the macroscopic A as algebraic functions of the relevant 
currents and also, for the latter, of the relevant vorticities as well as the electromag- 
netic field. Accurately realistic "equations of state" can only be obtained by detailed 
microscopic analysis, and even if available would be likely to be too complicated to be 
convenient in practice. Since use of a conservative model involving only three indepen- 
dent constituents is already a simplification, it is reasonable to employ a correspondingly 
simplified "equation of state" . 

A very convenient kind of simplification, that is implicitly or explicitly employed in 
most of the preceding work in a Newtonian framework || 0, and that will be adopted for 
the work of the present section, is to suppose that the action contribution (and hence the 
corresponding stress momentum energy contribution) of the "superfluid" or "quantised" 
part is separate from that of the "normal" or strictly "classical" part, meaning that at 
the level of the "mesoscopic" description the relevant "matter" contribution Am is taken 
to be decomposable as a sum in the form 

A M = A Q + A c , (79) 

where, in the three constituent models with which we are concerned here, the "quan- 
tum" part depends only on the neutron and the proton currents, while the "classical" 
part of the Lagrangian density depends only on the comoving electron and entropy 
currents, so that the variation (|7|) decomposes as the sum of two separate contributions 

8A Q = f p 5n? + < g v ° 5g pa , (80) 

and 

SA C = \i% K P - © Ss + < g v ° 8g pa . (81) 

It seems plausible that such a simplification should be fairly satisfactory for the low 
temperature limit, 0^0, thus providing an adequate approximation for many pur- 
poses in neutron star theory. However as was clearly understood by Landau when 
he obtained the equation of state for the "cool" limit of his original two constituent 
model |25[ (of which an elegant relativistic analogue has recently been provided by the 
present authors |26|]) there will inevitably be significant coupling between the superfluid 
and "normal" constituents at higher temperatures, for which the decomposition (|79|) 
would therefore cease to be valid. 

The simplest conceivable models would be based on not just partial but complete 
separation, which would mean that the contribution Aq in ( |79|) would itself split up 
as the sum of two distinct terms, A n and A p say, of which the first depends only 
on the magnitude of the neutron current vector, and the second only on that of the 
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proton current vector. However ever since the earliest pioneering work by Andreev 
and Bashkinpi it has been generally recognised that (even in the zero temperature 



limit) such an extreme simplification will be unrealistic. This means that the form of 
Aq will be such that the proton and neutron momenta will be given in terms of the 
corresponding currents by a (necessarily symmetric) matrix relation of the form 

lL* = gpa K r *nZ, (82) 

in which, as well as the diagonal terms /C nn and /C pp that would be present even in 
the oversimplified fully separated case, there will also be a cross term /C np = /C np . 
The presence of such a term implies that the respective directions of the neutron and 
proton momenta will deviate from that of the corresponding currents, in a manner that 
is familiar in the well known example of the Landau type two constituent model at 
non-zero temperature. Such a (non-dissipative) deviation effect is usually referred to as 
entrainment. (It is however to be noted that there is an alternative but regrettable usage 
whereby some authors call it "drag" , which is misleading, since the standard use of the 
latter term is for dissipative forces opposing relative movement.) The (symmetric) 2 by 
2 matrix relation ([82]) can be rewritten in inverted form as 

K = 9 P °K> Y ^% , (83) 

using the usual convention that (just as g pu denotes the components of the matrix 
inverse of the metric tensor g pa ) the matrix with (lowered) components K, T4> is the 
inverse of the original matrix with components /C r *. 

For the purpose of comparison with our treatment [^] of the single constituent case, 
for which it was found convenient to work with a "dilatonic amplitude", <& defined as 
the square root of the ratio of the particle number density to the effective mass, it is 
of interest to introduce analogous amplitudes <P n and <P P for the protons and neutrons 
respectively in such a way as to allow (R3T) to be rewritten more explicitly as 



rw = $1 ft + IC np (ft - ft) , rv = <P 2 p ft + IC Dp (ft a - ft a ) , (84) 
which evidently requires the specifications 

= ^nn + /C np , ^ = /C pp + /C np . (85) 

In order to relate the formula to (|82|) to the expressions traditionally used in the 
relevant literature as previously developed in a Newtonian framework [jTTj| , |28| , 0, it is 
convenient to introduce what from a relativistic point of view is a rather artificial 
concept, namely that of the "rest mass", m say per baryon, which can indifferently be 
taken to be the mass of a free neutron, the (more precisely definable) mass of a free 
proton, or as a compromise between the two, the mass of an ordinary free hydrogen 
atom in its ground state, or even the standard atomic unit as conventionally defined in 
terms of oxygen. The reason why it does not matter which of these one uses is that this 
"rest mass" is needed merely as a unit for calibration, and in any case by the standards 
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of accuracy of the present treatment all these various alternative definitions give are 
virtually the same result. However the calibration mass m is chosen, one can use it to 
define corresponding conserved "mass currents" p^ and p p for the protons and neutrons 
by the formula 

p v r = mn v r . (86) 

Another popular concept in the Newtonian context, despite the fact that its calibration 
is subject to a similar degree of arbitrariness, is that of what is commonly but mislead- 
ingly referred to as a superfluid "velocity". One can define such so called "velocities" 
v nu and v pu for the neutrons and the protons in terms of the chosen calibration mass 
m (as specified by one or other of the various conventions mentioned above for defining 
the "baryon rest mass") by setting 

Hi = mvi , (87) 

which means that v\ and v v v can be read out just as correspondingly rescaled momentum 
variables. It is to be observed that, no matter how the constant mass scale factor m 
is chosen, v nu and v vu can not be made to be "velocity" 4- vectors in the strict sense 
because their magnitudes will in general be variable, so that they can not be made 
to be unit vectors, except in such trivial special cases as that of a spatially uniform 
stationary state. In terms of these so called "velocities", which are actually just rescaled 
momentum variables, fl83|) can be rewritten in the form 

P^g^PvK, (88) 

where the coefficients p r4 , are components of a "density matrix" of the traditional but 
unsatisfactorily mass scale dependent kind, which will be given in terms of the more 
precisely defined matrix components K, r4> introduced above by 

P r0 = m 2 JC r4 , . (89) 

Writing this out using the definitions 

Pn = (m<P n ) 2 = p nn + p np , p p = (m<P p f = p pp + p np , (90) 

it can can be be seen that when translated into terms of the pseudo-velocities, the 
expansions (Q) will take the kind of form that is familiar in the Newtonian literature: 

Pna = PnVl+ Pnp« ~ <) , P P a = Pp V I + p np « ~ O ■ (91) 

6 Polarisation in partially separated models. 

For the purpose of extending the kind of "mesoscopic" model considered in the previous 
section to a "macroscopic" model representing the effect of averaging over a congruence 
of vortex lines in the manner described in Section ^, it is useful to start by using the 
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relations obtained in Section |5| to obtain rough estimates of the total electromagnetic 
flux, 

$ = j> A a dx u , (92) 

that is to be expected through a circuit at sufficiently large distance round a vortex 
defect characterised, according to (]51|) by neutron and proton and winding numbers 
u n and v v . It can be seen, by working out fl33| ) for the particular case of the proton 
current n° that the total value of the electromagnetic gauge field will be expressible in 
the form 

A a = -(< + «) + (93) 

e 

where the factor is a quantity that would be zero if there were no "entrainment" 
but that will in general be given by 

a£ = ^ = £=, (94) 
^pp Ppp 

while the final remainder term will be given by 
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Assuming that the coefficients /C r# tend approximately towards uniform values at large 
distance, the corresponding flux $ will tend approximately towards a limit given in 
terms of these values by 

$ ~ z/ r $ r + $ A , (96) 
where, by fl5T|) , the proton vortex flux contribution will be given simply by 

*p = ^, (97) 

which is just the usual flux quantisation unit associated with the relevant total charge 
2e of the Cooper type proton pairs, while the neutron vortex flux contribution would 
vanish if there were no "entrainment" but will in general be given by the formula 

$„ = —<, (98) 
e 



where is the relevant "entrainment factor". The formula (|98|), with given by 
the final expression in (@), is already familiar from the corresponding analysis in a 
Newtonian framework P7L EBJ. If the "Meissner effect" were fully effective the residual 



flux term would tend to zero, but as was originally realised by London in the context 
of ordinary metallic superconductivity this can be expected only in a strictly static 
background, whereas more generally even in a stationary case provided the background 
is rotating, there will be a residual magnetic field even at large distances from a vortex, 
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and hence a corresponding "London" flux contribution, which will be given in the 
present case by 



$ A = -—j>n£g p<7 dx c '. (99) 

Under reasonably stable conditions one would expect that energy would be locally 
minimised by avoidance of the field build up that would result from persistent devi- 
ation of the proton current contribution n° from the canceling background current 
contribution n° provided by the "normal" electrons, i.e. one would expect to have 

n° « n c u a , (100) 

where m m is the unit 4- velocity vector of the "normal" background, as introduced above. 
Under such circumstances, the "London" flux contribution ( j99|) can be expected to 
be given in terms of the asymptotic values of the proton number density n p or the 
corresponding conventionally normalised proton mass density as defined by 

u p = ( - n p V) 1/2 > Pp = mn P > ( 101 ) 
by an estimate of the form 



(J) £H i- 

A 



lu a dx u = lu a dx a . (102) 



However that may be, the "London" flux contribution $ A through a circuit round a 
plane surface element of area S can be attributed to an average orthogonally oriented 
magnetic field contribution (H) given by 

(H) = ^ . (103) 

In addition to this external London contribution, the corresponding total average mag- 
netic field, as defined by 

<*> = §, ( 104 ) 

will contain an internal contribution from the part of the flux that is confined within 
the neighbourhood of the vortex. According to (|96D, this internal contribution will be 
expressible by 

(B) - (K) = ±* r (v>T) , (105) 
Tin 

where the effective average neutron and proton vorticities (w n ) and (w p ) are given in 
terms of the corresponding quantised circulation integrals ( |5lD by 

(w Y ) = -jix r p dx p , (106) 

The formula ( |105| ) provides what one needs for going over from the "mesoscopic" 
level of analysis we have been using so far in this section to a "macroscopic" level of 
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analysis in terms of averages over many vortices. The average magnitude (w r ) is in- 
terpretable as the relevant component (as determined by the choice of orientation of 
the circuit under consideration) of the kind of macroscopic vorticity tensor that 
was considered in the preceding section. The total average magnitude (B) is analo- 
gously interpretable as the corresponding component of the macroscopic electromag- 
netic field tensor F p(T , while the external "London" contribution (H) will be similarly 
interpretable as the corresponding component of the generalised displacement field ten- 
sor H pa that was introduced in the previous section. By comparison with (|67D , it can 
thus be seen that the difference ( |105[ ) will be interpretable as the corresponding com- 
ponent of 47rAip a , where M. pa is the relevant polarisation tensor. It follows that this 
macroscopic polarisation tensor, representing the effect of the part of the magnetic field 
that is confined in the neighbourhood of the vortices, will be given in terms of the 
macroscopic proton and neutron vorticity tensors by 

An 2 hM pa = ^ r wj a . (107) 

It is important to point out that this result differs from the result obtained in the 
Newtonian limit by Mendell |7|], who found the polarisation to be zero. However that 
conclusion was based on the inclusion of an extra magnetic interaction term but on 
the neglect of a kinetic contribution that, according to a more recent analysis p9|| , will 
actually cancel it out. 

Working out the polarisation tensor explicitly in terms of the elementary proton and 
neutron vortex flux units ( p7[ ) and ( j98|) it can be seen that the Dirac-Planck constant 
Ti cancels out, leaving a relation of the form 

iireM pa = w% + alw a pa , (108) 

in which the last term on the right is entirely due to the "entrainment" effect, with- 
out which the neutron vorticity would not contribute to the magnetic polarisation. 
According to ( |57D the displacement field due to the "London effect" will be given by 

H^ = F p „-\w% + alw%). (109) 

In terms of the "normal" background's acceleration vector u v and rotation tensor 
Vt pv (whose magnitude Vt = (^f^f^") 1 / 2 is the local angular velocity) as defined by 

V [M u„] = 2VL pv - u [p u v] , ii M = u u V u , (110) 

one would expect that in steady circumstances such that the estimate ( |102| ) is valid, 
the "London" field would be close to a value given by 

TL TflD 

'Hp.y « -jp— (u|>tti/] - 2f2 AU/ ) = -{u [p u u] - 2Q pu ) , (111) 

in which it can be seen that the electric part is proportional to the acceleration while 
the magnetic part is proportional to the angular velocity. 
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The "London field" (|109|) will be interpretable (whether or not the conditions for 
( |1 1 1| ) are satisfied) as representing the contribution to the total macroscopic average 
electromagnetic field F pu from the magnetic field (if any) "outside" the vortices. Such 
a heuristic interpretation should not be taken too seriously because there will not be 
an absolutely clear cut boundary between the "inside" and the "outside" of a vortex. 
Nevertheless the distinction between the internal contribution described in terms of 
polarisation and the external contribution described as a "London" field will usually 
be fairly precise in practice, because the magnetic field can be expected to suffer an 
exponential decline with a lengthscale small compared to the intervortex separation. 

In so far as there is some residual degree of arbitrariness in the distinction between 
the internal "polarisation" contribution and the external "London" contribution, the 
formula ( |109| ) can be considered as an exact defining condition. However this convention 
will only be consistent with the original definition of the displacement field via (|65| ) in 
terms of partial differentiation with respect to the total field F pv if the Lagrangian 
density contribution A has the right functional form. For mathematical consistency 
between flB5|) and ( |109| ) it is necessary and sufficient that A should be decomposable in 
the form 

A = A MV + A F , (112) 

where the macroscopic contribution Amv is required to be functionally independent of 
F^, though it is dependent now on the vorticities w r as well as on the currents, while 
the remaining contribution is required to have the specific form 

k Y = ^-n pa n up ■ (ii3) 

Ion 

So long as the vortices are (as is expected to be the case in typical neutron star applica- 
tion) sufficiently far apart to justify the heuristic interpretation of Ti^ as the intervortex 
contribution to the macroscopically averaged field F pu , it is physically plausible that the 
action contribution of this "London" field should be given by substituting Ti prj for F^ in 
the standard Maxwellian formula ([12]). Since the formula ( |113| ) is precisely what results 
from such a substitution, the physical coherence of the description we have developed 
here is reassuringly confirmed. 

In the context of the partial separation approximation on which the analysis of this 
section is based, there will be no loss of generality in decomposing the macroscopic 
Lagrangian density contribution Amv i n the form 

A M v = A M + A v , (114) 

where Am is the function already used in the mesoscopic analysis, which in the present 
section is assumed to have the partially separated form (79), while Ay is an extra term 



which unlike Am is algebraically dependent not just on the currents but also on the 
vorticity 2-forms. 

Before concluding, it is to be remarked that, for the purpose of comparison with pre- 
vious work in Newtonian theory, it may be useful to introduce in our macroscopic model 
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the notions of "velocities" and of "density matrix" that were defined in a mesoscopic 
context in Section 5. This will mean to replace the (in our opinion more natural) forms 
([77]) and (|78| ) of the respective neutron and proton equations of motion by alternative 
decompositions 

2p a v»°V [rT <, = p np (v p ° - vn^K, + «CV P K p , (115) 

and 

2 Pp v^V [a vl } = p np (v™ - vn^-'K. ~ ^(Ppp^ + Pn P v na )F„v + <F P \; p , (116) 

in each of which the first term on the right can be interpreted as an "entrainment force" 
density. It is to be observed that - due to the antisymmetry of the generalised vorticity 
tensors- such entrainment forces always act orthogonally to the relevant relative differ- 
ence v ncr — v pcr , thereby respecting the strictly conservative character of the model. (The 
occasional use in the literature of the term "drag" instead of "entrainment" can be seen 
to be misleading, since a genuine drag force is intrinsically dissipative, and if present 
would not act orthogonally but on the contrary would be aligned with the relevant 
relative difference. In order to conform with standard aerodynamical terminology, the 
"entrainment" force should be described not as "drag" but rather as a kind of "lift".) 



7 Conclusions 

What remains for future work is the derivation of a suitable formula for the explicit 
form of the contribution Amv- Supposing that the vorticity independent mesoscopic 
contributions Aq and Ac are known, the "entrainment ratio" needed for the specification 
via ( |113|) of Ap will be given explicitly by @), and one might hope that (as in the 
example of our preceding study of the single constituent case|3[]) there might be a 
similarly simple and explicit formula for the corresponding form of the extra vorticity 
dependent term Ay. 

If the vortices were of purely "local" type, with not just their magnetic but also 
the associated circulating currents falling off exponentially in the exterior, then one 
would expect that Ay would just be a homogeneous linear function the magnitudes, 
w n = (^w^w^) 1 ^ 2 and w p = (^w^w^) 1 ^ 2 of the corresponding vorticity forms. In 
the absence of "entrainment" the elementary proton vortices would indeed be of this 
"local" type (like the electron vortices in ordinary "type II" metallic superconductors). 
However even in the single constituent case|| the neutron vortices are of "global" type, 
with energy logarithmically dependent on their separation, and in the present case the 



same will be true of the proton vortices, due to the non zero value that is expected p8 
for the entrainment factor a p given by (0). The simultaneous presence of neutron 
and proton vortices can give rise to interesting interaction effects, which have been 
examined in a Newtonian analysis by the Sedrakians||. However the most suitable 
way to represent these effects explicitly within the present framework is not yet clear 
to us. 
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Another important issue would be to work out the magnetohydrodynamic limit of 
our equations of motion in order to get rid of superfluous degrees of freedom and to 
obtain more tractable equations when one wishes to consider some specific problems 
such as glitches in neutron stars. 
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